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Neutron star models are studied in the context of scalar-tensor theories of gravity in the presence of a conformally coupled scalar 
field, using two different numerical equations of state (EoS) representing different degrees of stiffness. In both cases we obtain a 
complete solution by matching the interior numerical solution of the coupled Einstein-scalar field hydrostatic equations, with an exact 
metric on the surface of the star. These are then used to find the effect of the scalar field and its coupling to geometry, on the 
neutron star structure, particularly the maximum neutron star mass and radius. We show that in the presence of a conformally 
coupled scalar field, neutron stars are less dense and have smaller masses and radii than their counterparts in the minimally coupled 
case, and the effect increases with the magnitude of the scalar field at the center of the star. 
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